Exchange-correlation energy of a multicomponent two-dimensional electron gas 
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We discuss the exchange-correlation energy of a multicomponent (multi-valley) two-dimensional 
electron gas and show that an extension of the recent parametrisation 1J of the exchange-correlation 
energy by Attacalite et al. describes well also the multicomponent system. We suggest a simple mass 
dependence of the correlation energy and apply it to study the phase diagram of the multicomponent 
2D electron (or hole) gas. The results show that even a small mass difference of the components (e.g. 
heavy and light holes) decreases the concentration of the lighter components already at relatively 
high densities. 
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I. INTRODUCTION 

One of the most interesting discoveries in semicon- 
ductor heterostructures was the formation of a two- 
dimensional (2D) electron gas 2}. For example, this has 
set the starting point for building nanostructures such as 
quantum dots or wires (for reviews see 0,0, HI)- Many 
of the simplified models describing semiconductor nanos- 
tructures are based on the assumption of a homogeneous 
electron gas. The effects of the detailed band structure 
are included only by using an effective mass for the elec- 
trons and by screening the electron-electron interaction 
with a static dielectric constant. For the conduction elec- 
trons of a direct semiconductor, like GaAs for example, 
this approach describes fairly well the underlying elec- 
tronic structure 0,3- 

In the following, we call such a normal 2D electron 
gas as a two-component gas, the components being the 
spin-up and spin-down electrons. Similarly, a polarised 
electron gas is called a one-component Fermi gas. (In the 
literature, the latter is sometimes referred to as "spinless 
fermions" ) . 

Layered semiconductors, however, provide several 
structures where the simple picture of an ideal (one- or 
two-component) two-dimensional electron gas is likely to 
fail. In elemental semiconductors, silicon and germa- 
nium, the conduction band minimum consists of four 
equivalent minima at nonzero values of the A-vector 
and, moreover, these minima have non-isotropic effec- 
tive masses. The resulting conduction electron gas has 
an "internal" degeneracy of eight (two from spin and four 
from the four minima), and is a possible example of an 
eight-component Fermi gas. The correlations in such an 
electron gas are much more complicated and have been 
studied in connection with the 3D electron-hole plasma 
in semiconductors [||. 

The valence band maximum in most semiconductors 
is degenerate, consisting of two bands corresponding to 
heavy and light holes. In 3D structures the hole gas will 
then consist of two kinds of holes. In 2D structures there 
will be an energy separation between the band minima 
due to the fact that heavy holes can be more easily lo- 



calised than the light hole. Nevertheless, if the mass dif- 
ference is not too large, both kinds of holes can coexist 
and the correlations between them have to be consid- 
ered. The hole gas of most semiconductors is one exam- 
ple of a four-component Fermi gas. Electron layers in 
semiconductor heterostructures provide yet another sys- 
tem of multicomponent electron (or hole) gas, intensively 
studied theoretically 0, H, 0, El- Recently, electron ad- 
dition spectra of vertical quantum dots based on layered 
structures have also been measured [Tl]. 

Density functional theory (DFT) in the local (spin) 
density approximation (LSDA) has shown to be a use- 
ful starting point for studies of the ground state proper- 
ties of electrons in quantum dots0. The 2D conduction 
electron gas of the semiconductor is then treated as an 
interacting electron gas. In the LSDA, the exchange- 
correlation energy is approximated by that of the ideal 
2D electron gas, the effective mass and the static dielec- 
tric constant appearing only as scaling factors (for energy 
and distance). Applying the density functional method 
to the above mentioned multicomponent structures calls 
for a parametrisation of the exchange-correlation en- 
ergy (e xc ) for the multicomponent system in question. 
Tanatar and Ceperley0 computed the one and two- 
component (polarised and non-polarised 2D electron gas) 
using the quantum Monte Carlo technique. More re- 
cently Attaccalitc et presented new Monte Carlo 
results for these systems. They also derived an inter- 
polated expression for the exchange-correlation energy 
functional to be used in DFT calculations. The most im- 
portant difference to previous parametrisations like for 
example the one by Tanatar and Ceperly^^l is that here 
the Monte Carlo-calculations are performed also at inter- 
mediate polarisations. The usual approximation to treat 
the polarisation dependence of the correlation energy on 
the same footing as the (exactly known) exchange part 3 
can thus be abandoned. 

For a multicomponent system the situation becomes 
more complicated and the transferability of e xc from one 
system to another is not so obvious. The purpose of 
this paper is to estimate the exchange-correlation en- 
ergy of a general multicomponent 2D fermi gas. We will 
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first discuss the general high and low density limits and 
the effect of mass difference between the components. 
We will then show that an analytic continuation of the 
parametrisation by Attaccalite et aL0 fits well the exist- 
ing many-body calculations for multicomponent systems. 
As reference data we use the calculations of Conti and 
Senatore^(j for a four-component 2D electron gas and 
the results of Apaja et a/.|l3| for a 2D charged Bose gas 
which can be viewed as an upper bound for the energy of 
a Fermi gas with an infinite number of internal degrees 
of freedom. 

Finally, we will study the phase diagram suggested for 
this energy functional and show that already a small mass 
difference will decrease the concentration of the lighter- 
mass components. When the density gets low enough, 
eventually all systems become one-component, as pre- 
dicted by the Attaccalite et al. parametrisation for the 
normal 2D electron gas. 



II. EXCHANGE-CORRELATION ENERGY: 
GENERAL DISCUSSION 

We consider a homogenous gas consisting of A different 
kinds of fermions, or components, as we call them. For 
instance, the different spin states of the homogenous gas 
are treated as different components. For a normal elec- 
tron gas A = 2 and for a fully polarised gas A = 1. We 
assume that all the components have the same charge, 
but they may have different masses. The total density n 
of the gas is 

A A 

n = y"Vj = ny^i/j, (1) 

i=l i=l 

where is the number density and vi = rii In is a dimen- 
sionless concentration of the component i. Note that the 
concentrations are normalised as Vj = 1. The den- 
sity parameter r s — 1/ 1 \fwn is always related to the total 
density n of the gas. The kinetic energy per particle of 
the noninteracting gas becomes 

e^W,W) = ?E^ (2) 
ri * — ' to,- 

where TOj is the mass of the component i. The total 
exchange energy is the sum of exchange energies of the 
different components. The total exchange energy per par- 
ticle is then 

e x (r s ,{^})--^—^J2^ /2 - ( 3 ) 

Note that while the kinetic energy depends on the masses 
of the components, the exchange energy is independent 
of the masses. The correlation energy per particle is de- 
fined as the difference between the exact energy and the 



Hartree-Fock energy, e c = e to t — (efe + e^)- In true semi- 
conductor systems the different components can repre- 
sent, for example, different band minima. The exchange 
energy has then also components coming from the ex- 
change between electrons belonging to these different 
minima. Generally such exchange interaction is found 
to be small [l5l[l^|. In the spirit of the effective mass ap- 
proximation we will completely neglect such effects here. 

In order to study the multicomponent correlation en- 
ergy we will first consider the case where all the compo- 
nents have the same mass. In this case the exchange cor- 
relation energy can be written as a function of the density 
parameter r s and the concentrations vi. The symmetry 
requires that the energy functional is symmetric with re- 
spect of the interchange of the different concentrations 
Vi, i.e., it does not depend on the indexing. We define 
numbers 

^y = X>7, (4) 

i=l 

which by construction have the desired symmetry. We 
notice that Z\ = 1, e x oc Z3/2 and the kinetic energy ek <x 
Zi (for equal masses). In general any function g{{v{\) 
with the desired symmetry property can be written as 
a function of the ^-numbers, for example, as a series 
expansion 

g{{ p i\) = f{{ z ~a}) = aiZ 7l + a 2 Z l2 + a i2 Z 7l Z 72 H , 

(5) 

where are constants and 7^ can be noninteger. One 
way to approach the correlation energy is to use the 
suggestion of Wigner|l4| who interpolated the exchange- 
correlation energy of an electron gas between the two lim- 
iting cases, the pure exchange in the high density limit 
and the energy of the Wigner crystal in the low density 
limit: 

£xc(r s , {vi}, {rrii\) = e W c(r a ) 

(6) 

+g(r s , {>;}, {m l })[e x (r s , {v.j}) - e wc (r s )] 

where ewe is the classical Madelung energy of the 
Wigner crystal, which does not depend on the masses. 
The advantage of this approach is that the mass de- 
pendence appears only in the interpolation function g. 
Furthermore, if the masses are equal, the interpolation 
function g can immediately be written as in Eq. JSJ. 

In order to estimate the interpolating function g one 
would need results of accurate many-body calculations 
for systems with different numbers of components, and 
for several values of r s and the concentrations . We are 
not in the position to perform such computations at this 
point but use existing data for special cases. For the two- 
component system Attaccalite et a/.Q have performed 
Monte Carlo calculations. They fitted an accurate inter- 
polation formula to the results which is valid for any r s 
and polarisation. 

For the four-component system, Conti and 
Senatore^(j have presented Monte Carlo results 
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for several values of r s in the case where all the 
concentrations are equal. For systems with more 
components there exist no accurate data. However, the 
limit of infinite components can be approximated by 
the results obtained for charged bosons |l3j. If all the 
concentrations are the same (yi = 1/A) in the limit 
of infinitely many components, all the particles are at 
different internal state and the Pauli exclusion principle 
can not prevent putting all the particles in the same 
orbital state. The symmetric boson state is then a 
legitime state for such a fermion system and provides 
an upper bound for the energy. However, the existing 
results beyond the one and two-component systems are 
still quite limited for an accurate interpolation of the 
function g. 

A. Mass dependence 

As discussed above, the exchange energy is mass- 
independent and even in the multi-component electron 
gas only the kinetic and the correlation energies will de- 
pend on the masses. If all the components have the same 
mass M, the kinetic energy of Eq. @ can be written as 
€k = h 2 Z 2 /(r 2 M). In the general case, the kinetic en- 
ergy can still be written in the same form by defining an 
average mass M as 

1 i A 2 

- = -y^. (7) 

M Z 2 ^ ra, V ; 

We notice a scaling: By deviding the the kinetic and 
exchange energies by M, the resulting ratios will only 
depend on the product Mr s (e k /M = h 2 Z 2 /{Mr s ) 2 and 
e x /M = — 8e 2 Z 3 / 2 /(47r£oe37rMr s )). Moreover, consider- 
ing the atomic units in the case that all masses are the 
same, we notice that (in this case) also the total energy, 
and thus the correlation energy has the same scaling. It 
is then a reasonable first approximation to take the mass 
as an average mass and use it as a scaling factor 

M 

£xc{r s , {vi\, {mi}) = — e IC (¥r s ,{i/,}, {mi = m e }) , 
m e 

where m e is the bare mass of the electron (or any suitably 
chosen effective mass). (Note that within this approxi- 
mation, the interpolation function g in Eq. © would be 
written as g(Mr s , {Z^})). 

III. EXTENSION OF THE TWO-COMPONENT 
FUNCTION 

For estimating an interpolation formula for the mul- 
ticomponent correlation energy, we start by considering 
the parametrised form of the two-component 2D elec- 
tron gas by Attaccalite et al. Q], in order to see if their 
function can be written in terms of the r s and the Z 7 - 
numbers, as discussed above. Attaccalite et al. use the 



common description where the energy is written as a 
function of r s and the polarisation £ = (nj — n 2 ) /n — 
v\ — v 2 of the electron gas. They chose a correlation 
energy function e c (r s ,C) = e xc (r s ,C) - £x(r s ,C) which 
diminishes the contributions of e x beyond fourth or- 
der in £ as r s increases. The parametrisation satisfies 
known high- and low-density limits: For high densities 
i.e. small r s the exchange-correlation energy obeys the 
perturbation theory result e xc = e x + ao(C) + bo(C)r s hi?V 
For low densities the Wigner crystal limit is recovered, 
e xc ~ ~ci/r s + c 2 /r^ 2 , where the constants c\ and c 2 
are independent of the number of components. 

Using the concentrations Vi, the square of the po- 
larisation £ of the two-component gas becomes C 2 — 
2(y\ + v 2 ) — 1 = 2Z 2 — 1. Now the exchange-correlation 
energy introduced by Attacalite et al. can be written as 
(using their notations except for the above replacement 
for C) 

e xc (r s , {vi}) = e~P r ° [e x - el 6) ] 

(8 

+e ( x 6) + a (r s ) + «i(r s )(2Z 2 - 1) + a 2 (r s )(2Z 2 - l) 2 , 

where e ( x 6} = (1 + f (2Z 2 - 1) + ^ (2Z 2 - l) 2 )e x (r s , C = 0). 
The functions tti(r s ) are parametrised by Attacalite et 
aL[l| and are independent of the concentrations. This 
two-component exchange-correlation function of Attac- 
calite et al is thus already of the desired form: It depends 
only on r s , Z 2 and Z 3 / 2 (through e x ) and can be directly 
extended to any number of components. 

For a given value of r s Attaccalite et aZ0] fitted their 
function to several values of polarisations, Q S [0,1], 
which correspond to Z 2 € [0.5,1] and Z 3 / 2 £ [v2>l]- 
An analytic continuation to larger number of components 
means an extrapolation of the e xc to the region Z 2 € [0, 1] 
and Z^j 2 £ [0, 1] for a system with infinitely many com- 
ponents, but only to Z 2 £ [0.25, 1] and Z 3 / 2 € [0.5, 1] for 
a four-component system. Moreover, since Z 3 / 2 only ap- 
pears in the exchange energy, Eq. @, which is exact for 
any number of components, errors in the extrapolation 
only result from the extension of the Z 2 space. 

The multicomponent gas must obey the low density 
Wigner crystal limit and the high density exchange limit; 
both of these requirements are satisfied by Eq. @. By 
taking a "paramagnetic" gas, i.e. by setting Vi = A -1 for 
all the components and fixing r s , it can be seen that the 
exchange-correlation energy above is also a monotonous 
function of the number of components. 

Figure 1 shows the exchange-correlation energy derived 
from Eq. @ for selected values of r s as a function of Z 2 . 
For the interval Z 2 £ [.5, 1] the lines correspond to the 
normal two-component electron gas with partial polarisa- 
tion. For the interval Z 2 £ [0, .5] the lines are derived for 
the case where all the concentrations are equal, Vi = 1/A 
and consequently Z 2 = 1/A. The points are results of 
different many-body calculations. Note that the lines 
from Z 2 = 0.5 to Z 2 = 1 agree exactly to the results of 
Attaccalite et alW. We can see that the extension to a 
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FIG. 1: Dependence of the exchange-correlation energy of the 
multicomponent electron gas as a function of parameter z?2, 
defined by Eq. (JIJ. The lines show the result of the inter- 
polation formula of Attaccalite at aZ.Q, Eq. (|5J, extended 
to multicomponent systems. The points correspond to nu- 
merical results of many-body calculations as follows: Black 
dots for Z2 — 1 and Z2 = 0.5 are from Ref.Q; black dots for 
Z 2 = 0.25 from Ref.[l3hblack dots for Z 2 = from Ref.0; 
open circles from Ref.|l2j|: black squares from Ref . [TtI . Re- 
sults for r s — 2, 10, and 20 are shown. 

four component system (Z2 = -25) is fairly accurate and 
the extension even to the boson case {Z 2 — 0) is still 
reasonable. In fact, the agreement of the Attaccalite el 
aL-function with the existing many-body results for four- 
and infinite-component systems is so good that there is 
no need to make a better interpolation formula for the 
multicomponent system, before more total energy calcu- 
lations are available in the region Z2 < .5. It would be 
important to have results also for cases where the con- 
centrations of all the components are not equal. 



TABLE I: Estimated values of r s where a A-component sys- 
tem spontaneously transforms to a one-component system, 
assuming that all the components have the same mass. 
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FIG. 2: Phase diagram for the four-component gas in the 
rs-rrih/mi plane. The solid lines separate the 1, 2 and 
4-component phases (denoted bu numbers). In the four- 
component phase the concentrations of the components with 
the heavy mass, say V3 and V4, decrease continuosly from 0.5 
to 0.25 in going from the lowest solid line to the mass ratio 1 
or to r s — 0. The dashed lines correspond to concentrations 
vz = v& = 0.42 and 1/3 = V4 = 0.35. 



IV. PHASE DIAGRAM OF THE 
MULTICOMPONENT 2D ELECTRON GAS 

It was recently suggested by Attacalite et al.\j^ that 
the 2D electrongas shows a transition to a polarised 
gas at r s « 25 We will first study the stability of 
the multicomponent gas assuming all the masses to be 
equal. The extended exchange-correlation energy func- 
tional suggests that all multicomponent cases have a 
transition to a one-component phase when r s increases. 
Table I gives the estimated transition points. For any 
number of components, the transition happens nearly at 
the same point and directly from the multicomponent to 
a one-component phase. It is caused by the slightly dif- 
ferent r s -dependence of the one-component total energy 
as compared to the others. For A > 2 the total energy is 
a monotonously decreasing function of A (for any r s ). 

Next we will study the effect of the mass difference, 
using the approximation of Eq. © . We assume a four- 
component system and fix the masses pairwise equal so 



that two components are heavier than the two others, 
say mi = rri2 = mi and to 3 = m 4 = rrih- Note that due 
to the spin degeneracy of electrons (or holes) there will 
always be two components with the same mass, and the 
number of components will be even, except of the spe- 
cial case one. The phase diagram of the four-component 
electron gas obtained using Eq. © is shown in figure 
2. We show the maximum of the concentrations (1/4) 
as a function of r s and the mass ratio m = rrih/mi- 
When the masses are equal, also the concentrations re- 
main equal {i>i = 0.25) until r s ps 26.2 at which point 
the system becomes a one-component gas. When the 
mass ratio increases, the concentrations ^3 and V4 in- 
crease due to the reduction of the kinetic energy, and 
eventually vj, = 1/4 = 0.5 and the system has become 
two-component. A further increase of the mass ratio 
eventually changes the system to a one-component gas. 
We notice that as the mass is increased the polarisation 
occurs at smaller r s . The reason is the scaling of the total 
energy and the dimensions with the mass. 

In real semiconductor heterostructures the mass differ- 
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ence can arise from different constituents of the two layers 
in a double layer system or, in the case of holes, simply 
from the mass difference of the heavy and light hole. In 
both cases, however, there can be a constant potential en- 
ergy difference from one component to another. This en- 
ergy difference arises from the localisation of the particle 
in the direction perpendicular to the 2D-layer. For exam- 
ple, if we consider holes in a ID harmonic potential, the 
ground state energy (the lowest perpendicular mode) has 
the energy eo(m) = Hloq/ 1 {2^Jm/m e ), where ujq is the con- 
finement strength for m e . This energy is smaller for the 
heavy hole and, consequently, further favours the transi- 
tion to a two-component system where only heavy holes 
exist. In order to have a 2D electron gas the Fermi energy 
has to be clearly below the first excited perpendicular 
state, ei(m) = Zhu>o/{2y/m/m e ) in a harmonic well. In 
order for the heavy (mass m^) and light holes (mass m{) 
to coexist, the ratio (e^mi) ~ e (m h )) / (ei(m h ) - e (m h )) 
has to be clearly smaller than one. For example in the 
cases of Si, Ge and GaAs (m^/m/ > 5) this ratio is larger 
than 0.6. Consequently, a 2D layer with both heavy and 
light holes is not possible. However, there are semicon- 
ductors, like AlAs, where m/i/mj ~ 2 and the above en- 
ergy ratio is only about 0.2, making the four-component 
gas possible. 



be used in density functional calculations. After general 
considerations of the concentration and mass dependence 
of the correlation energy, we have demonstrated that the 
recent parametrisation by Attacalite et al.\j^ of the two- 
component exchange-correlation energy has functional 
properties which allow a direct extension to multicom- 
ponent systems. Furthermore, although the function is 
fitted only to the two-component electron gas, it gives 
fairly accurate results for four-component and infinite- 
component systems. 



We suggest that the mass dependence of the corre- 
lation energy can be estimated with a properly cho- 
sen average mass. The application for multicomponent 
2D gas shows an interesting phase diagram. For equal 
masses all systems transform from a multicomponent di- 
rectly to a one-component gas when r s increases to about 
25 • • • 27 depending on the number of components. Even 
a small mass difference decreases the concentration of the 
lighter component and it transforms first to a purely two- 
component and then to a one-component system when r s 
increases. 



V. CONCLUSIONS 

A more accurate interpolation formula requires exten- 
We have studied the possibility to formulate a gen- sive many-body calculations for multi-component gases, 
eral multicomponent 2D exchange-correlation energy to We hope that these considerations encourage such work. 



[1] C. Attaccalite, S. Moroni, P. Gori-Giorgi and G.B. 
Bachelet, Phys. Rev. Lett., 88, 256601 (2002). 

[2] T. Ando, A.B. Fowler, and F. Stern, Rev. Mod. Phys. 
54, 437 (1982). 

[3] T. Chakraborty, Quantum Dots: A survey of the prop- 
erties of artificial atoms (North-Holland, Amsterdam 
1999). 

[4] S.M. Reimann and M. Manninen, Rev. Mod. Phys. 74, 
1283 (2002). 

[5] W. G. van der Wiel, S. De Franceschi, J. M. Elzerman, 

T. Fujisawa, S. Tarucha, and L. P. Kouwenhoven, Rev. 

Mod. Phys. 75, 1 (2003). 
[6] P. Vashista, P. Bhattaracharyya and K. S. Singwi, Phys. 

Rev. B 10, 5108 (1974). 
[7] D. Neilson, L. Swierkovski, J. Szymariski, and L. Liu, 

Phys. Rev. Lett. 71, 4035 (1993). 
[8] S. Das Sarma and P.I. Tamborenea, Phys. Rev. Lett. 73, 



1971 (1994). 

[9] M. Alatalo, P. Pietilainen, T. Chakraborty, and M. 
Salmi, Phys. Rev. B 49, 8277 (1994). 
[10] S. Conti and G. Senatore, Europhys. Lett., 36, 695 
(1996). 

[11] D.G Austing, T. Honda, K. Muraki, Y. Tokura, and S. 

Tarucha, Physica B 249-251, 206 (1998). 
[12] B. Tanatar and D.M. Ceperley, Phys. Rev. B, 39, 5005 

(1989). 

[13] V. Apaja, J. Halinen, V. Halonen, E. Krotscheck, and M. 

Saarela, Phys. Rev. B, 55, 12925 (1997). 
[14] E. Wigner, Phys. Rev. 46, 1002 (1934). 
[15] D. Beliaev, L.M. Scolfaro, J. R. Leite, and GM. Sipahi, 

Phys. Stat. Sol. (b) 2 10, 777 (1998). 
[16] Y. Hada and M. Eto, |cond-m at /0304228 (2003). 
[17] A. Gold, Z. Phys. B 89, 1 (1992) 



